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ABSTRACT 

The luminosity function for contact binary stars of the W UMa-type is evaluated 
on the basis of the ASAS photometric project covering all stars south of S = +28° 
within a magnitude range 8 < V < 13. Lack of colour indices enforced a limitation 
to 3373 systems with P < 0.562 days (i.e. 73% of all systems with P < 1 day) where 
a simplified My (log P) calibration could be used. The spatial density relative to the 
main sequence FGK stars of 0.2%, as established previously from the Hipparcos sample 
to V = 7.5, is confirmed. While the numbers of contact binaries in the ASAS survey 
are large and thus the statistical uncertainties small, derivation of the luminosity 
function required a correction for missed systems with small amplitudes and with 
orbital periods longer than 0.562 days; the correction, by a factor of 3x, carries an 
uncertainty of about 30%. 
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1 INTRODUCTION 



ASAS - All Sky Automated Survey 1 , is a long term 
project dedicated to detection and monitor ing variability of 
bright stars using small telescopes. Recentlv lPaczvriski et all 
( 2006) presented results of the V-filter observations obtained 
at the Las Campanas Observatory with a single 7 cm tele- 
scope. More than 50 thousand variables brighter than 14 
magnitude were discovered in about three fourth of the sky 
with S < +28°, among them 5,384 contact binary stars. 

This paper addresses the matter of the spatial density 
of contact binaries of the W UMa type (P < 1 day, from 
now on called "EW") as indicated by the ASAS survey data 
within its magnitude limits of 8 < V < 13. The relative spa- 
tial density (in relation to the main sequence FGK dwarfs) is 
evaluated via the luminosity function for the absolute mag- 
nitude range 1.5 < My < 5.5. 

We are very far from having a volume limited sample 
of EW binaries. In fact, the main problem is the complete- 
ness to a given magnitude limit. At present, the only sample 
which appears to include all EW binaries with photometric 
variability amplitud es > 0.05 mag , is the shallow Hipparcos 
sample to V = 7.5 jRucinskill20021) (from now on called the 
"Hipparcos sample" or the "7.5 mag. sample") consisting of 
35 objects, with 1/3 in that number being new detections. 
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1 See: [http://www.astrouw.edu.pl/~gp/asas/ asas. html and 
|http://archive. princeton.edu/~asas/| 



This sample suggests the relative spatial density of the EW 
binaries at a level of 0.2% of the FGK main sequence stars, 
but this fraction requires a confirmation by a larger sample, 
particularly for My > +3.5 where the Hipparcos sample 
contained very few objects. 



2 AMPLITUDE DISTRIBUTION AND 
MISSING SYSTEMS 

Spatial density estimates of variable stars are crucially de- 
pendent on the detection limitations: Photometric variabil- 
ity must be some 5 times larger than the photometric error, 
a, for detection while another margin of > 10 a is needed 
for proper characterisation and assurance that the star is 
not one of rather common S Set or Dor short-period pul- 
sating m ain-sequence st ars <Rucinski2 004l. Theoretical pre- 
dictions jRucinskl200lD indicate that low-amplitude contact 
binaries should be very common as the amplitude distribu- 
tion is expected to rise for amplitudes going to zero. 

EW binaries discovered in the ASAS survey show a dif- 
ferent amplitude distribution from that of the Hipparcos 
sample. Figure gives the distribution for the Hipparcos 
sample (left panel) and for the ASAS sample (right panel), 
both compared with the theoretical amplitude distribution 
expected for random orbital incli nations and for a flat mass- 
ratio distribution, Q(q) = const iRucinskill200lT) . The latter 
assumption is made for lack of any evidence to the contrary 
and needs verification. Also, both observational samples in- 
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Figure 1. The amplitude distributions for 35 contact binaries 
of the Hipparcos 7.5 mag. survey (the left panel) and for 4640 
binaries of the ASAS survey (the right panel) with P < 1 day 
are shown together with the expected distribution for randomly 
distributed orbital inclinations and a flat distribution of mass 
ratios (the dashed line). In both panels the theoretical curve is 
normalized to the total number of systems, although we show it 
also scaled up by 2 X and 3x in the right panel. 

elude a small admixture of contact binaries with unequally 
deep eclipses (sometimes called EB); the presence of these 
binaries tends to elevate the high amplitude end of the dis- 
tribution, so that the agreement of the Hipparcos and the 
theoretical distributions may be fortuitous. 

It is easily noticeable in Figure that while the Hip- 
parcos sample follows the theoretical prediction, but suffers 
from the low number statistics., the ASAS distribution lacks 
systems with amplitudes < 0.4 mag. This can be understood 
in terms of the detection selection effect with a relatively 
high detection threshold. For a direct comparison of the Hip- 
parcos and ASAS samples, the numbers estimated from the 
latter must be multiplied by a correction factor of about 
2 — 2.5 x (see Figure 0, with a considerable uncertainty of 
some 25% or so. 

Compared with the theoretical and Hipparcos distri- 
butions, the ASAS amplitude distribution appears to have 
some excess of systems with amplitudes in the range of 
0.4 < a < 0.6 and too few systems with amplitudes > 0.6 
mag. The absence of large amplitudes may be in the pho- 
tometric blending of images and thus a "dilution" of the 
amplitudes because of the large images given by the 7 cm 
telescope, particularly when coupled with the now recog- 
nized prevalence (> 50%) of triple and multip le systems 
with contact binaries iPribulla fc RucinskilfeOOa) . 



Table 1. The sample used to establish My = My (log P calibra- 
tion. VW Cep is the only system where corrections for the third 
component with L3/L12 = 0.08 has been applied. The last four 
systems were not used, see the text. 



Name 




B - V 


E B -v 


P(d) 


M v 


eM v 


VW Cop 


7.38 


0.82 


0.00 


0.2783 


5.17 


0.06 


OU Ser 


8.25 


0.62 


0.01 


0.2968 


4.41 


0.13 


SX Crv 


8.95 


0.52 


0.04 


0.3166 


4.01 


0.24 


YY Eri 


8.16 


0.66 


0.00 


0.3215 


4.43 


0.14 


W UMa 


7.76 


0.62 


0.00 


0.3336 


4.28 


0.11 


GM Dra 


8.65 


0.48 


0.03 


0.3387 


3.59 


0.19 


V757 Ccn 


8.40 


0.65 


0.02 


0.3432 


4.10 


0.17 


V781 Tau 


8.55 


0.58 


0.02 


0.3449 


3.94 


0.24 


GR Vir 


7.81 


0.56 


0.01 


0.3470 


4.15 


0.14 


AE Phe 


7.56 


0.64 


0.00 


0.3624 


4.12 


0.09 


YY CrB 


8.64 


0.62 


0.01 


0.3766 


3.89 


0.16 


V759 Con 


7.47 


0.53 


0.02 


0.3940 


3.41 


0.13 


EX Leo 


8.17 


0.53 


0.00 


0.4086 


3.13 


0.24 


VS66 Oph 


7.47 


0.41 


0.03 


0.4096 


3.10 


0.17 


AW UMa 


6.84 


0.33 


0.01 


0.4387 


2.71 


0.13 


CN Hyi 


6.57 


0.41 


0.01 


0.4561 


2.72 


0.08 


TY Men 


8.11 


0.29 


0.01 


0.4617 


1.94 


0.23 


RR Cen 


7.32 


0.34 


0.03 


0.6057 


2.17 


0.19 


IS CMa 


6.87 


0.31 


0.01 


0.6170 


1.84 


0.15 


V535 Ara 


7.15 


0.30 


0.02 


0.6293 


1.83 


0.22 


S Ant 


6.29 


0.30 


0.01 


0.6484 


1.88 


0.12 
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Figure 2. The period - luminosity calibration for the 21 EW 
systems with good Hipparcos parallaxes and free of triple-system 
complications. The open and closed circles mark systems bluer 
and redder than (B — V)o = 0.5, respectively. The line is given 
by: My = —1.5 — 121ogP with the orbital period P in days. It 
was derived from 17 systems with logP < —0.25. 



3 ABSOLUTE MAGNITUDE CALIBRATION 

The previously p roposed absolute magnitude calibration for 
contact binaries jRucinskil994l : lRucinski fc Duerbeckll997t 
iRucinskilliooi) . is a linear one: My — ap log P + o,bv(B — 
V)o + do, where (B — V)o is the reddening corrected colour 
index (which can be replaced by say (V — I)o); it can predict 
My to about 0.2 - 0.25 mag. It has been used successfully for 
many applications, particularly for confirming or d isproving: 
membershi p of contact binaries in stellar clusters jRucinskil 
1998, 2000). T he main basis for the calibration is the Hip - 
parcos sample tRucinskilll994l : lllucinski fc Duerbecklll997|) . 
This sample has been carefully scrutinized for multiplicity 
and for small, but measurable amounts of the local inter- 
stellar (IS) reddening. The former requirement is related to 



the high frequency of triple systems among contact binaries 
while the latter stems from the availability of sensitive indi- 
cators of small amounts of local reddening such as the linear 
polarization, EUV extinction or IS cirrus emission. 

The calibrating sample is given in Table Q It consists 
of 21 EW from the Hipparcos 7.5 mag. sample llRucinskil 
2002). Except for for the orbital periods (known practically 
without any error) and for the parallax and parallax error 
data taken from the Hipparcos mission catalogue (not listed 
here, but available in previous tabulations), the remaining 
data have been modified: The known triple systems have 
been removed (with the exception of VW Cep, which is im- 
portant for low My; its ratio Z/3/Z/12 is well known) and the 
magnitudes and colour indices were corrected for newly de- 
rived values of the IS reddening. The reddening is given to 
0.01 mag. for consistency with typical accuracy of B — V, but 
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its value is frequently known much better. The table should 
be regarded as an improvement over the previous sim ilar 
tabulations (|Rucinski fc DuerbedJll997t iRucinsriEoO^l in 
the treatment of the IS reddening and in the careful screen- 
ing for triple components. 

The ASAS data presented bv IPaczvhski et alJ ll2006l) 
are in the V band. The colour indices are currently un- 
known, but will be available in the future. We derive here a 
simplified My = My (log P) calibration utilizing the known 
loose correlation of th e periods and colour indices, as no- 
ticed long time ago bv lEggenl <ll967l) . We claim that a very 
simple relation of the form My = ao + ap log P may suf- 
fice in situations like the current one when an additional 
limit is added on the length of the orbital period. The very 
steep dependence is indeed visible in Figure Some scatter 
is due to a range in the reddening corrected colour index, 
0.28 < (B — V)o < 0.82; the spread is noticeable in the fig- 
ure where different symbols are used for (B — V)o above and 
below 0.5. Four binaries with logP > —0.25 are all blue, so 
a decision was made to use only 17 binaries with periods 
shorter than this limit. By using the bootstrap sampling 
technique of repeated least-squares solutions, we derived a 
simplified relation: My = -1.5(±0.8) - 12.0 (±2.0) log P, 
with a — 0.29 for a single object. The errors of the coef- 
ficients have been estimated from the bootstrap sampling 
distributions encompassing 67% of the results. 

If a simplified calibration works, a question arises 
whether the colour term is needed at all. The data have 
well established errors eMy derived from the parallax errors, 
so one can check if x 2 = X] — C) 2 /a 2 is approximately 
equal to the number of degrees of freedom. The 3-term lin- 
ear fit My = ao + ap log P + ap v (B — V) gives Xs\ = 13-4 for 
v = m — 3 = 14 degrees of freedom, thus indicating a basi- 
cally perfect fit and the need of a 3-term fit 2 . The simplified 
2-term linear fit, My = -1.5 - 12.0 log P, gives xl = 52.9 
for v = m — 2 = 15 degrees of freedom, confirming that 
indeed three terms are needed. However, we can invert the 
problem and ask if the decrease in x 2 due to the addition of 
the colour term is significant. It appears that the difference 
betw een x 2 and x 2 , tested with the F Snedeckor distribu- 
tion taamiltonll964l) suggests a probability of an accidental 
improvement at the level of 37%, quite a distance from the 
normally desired < 5%, so it would appear that the third 
term is not needed after all. Thus, the situation is somewhat 
confused. However, a consideration of the involved physics 
strongly indicates that a colour term should be present in 
the calibration unless it is not available, as in the current 
application. 



4 DISTANCES AND SPATIAL DENSITY 

The limitation to periods log P < —0.25 or P < 0.562 days 
means that from now on we will consider only the EW bina- 
ries from the main peak of period distribution, as shown in 

2 The coefficients for this sample are: ap = — 5.8 ± 1.1, apy = 
+2.75 ± 0.45, ao = —0.31 ± 0.32. Because of the period — colour 
(ap — apy) correlation, the My predict i ons a re practically the 
same as for coefficients in iRucinskil iiool l2004h . For consistency 
with the previous results, it is recommended to continue using 
the old coefficients. 
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Figure 3. The period distribution for the ASAS sample of 4640 
EW systems with P < 1 day. The My calibration derived here 
applies to logP < —0.25 and this limit is marked by the vertical 
dotted line. 
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Figure 4. Distances to 3377 EW systems with logP < -0.25 
plotted versus My. The distance limits for V = 8 and V = 13 
are shown by lines. 

Figure |3 Thus, out of the full sample of 5384 EW binaries 
in the ASAS survey, from now on we will discuss the results 
for 3377 systems or 63% of the whole sample or 73% of the 
sub-sample of 4640 systems with P < 1 day. For a direct 
comparison of numbers with the Hipparcos sample, a cor- 
rection factor of 1/0.73=1.37 must be applied. This factor, 
combined with the under-counting due to loss of low am- 
plitude binaries (Section |5J implies that the ASAS counts 
must be multiplied by a factor of about 3 x carrying an un- 
certainty of about 30%. 

With the known values of My, the distances are simply 
calculated from d = dex((V — My + 5)/5), where we ne- 
glect the interstellar reddening. As we can see in Figure [I] 
the distances are moderate, so this assumption is accept- 
able except for the bins of higher luminosities where the 
distances extend beyond a few hundred parsecs. Thus, we 
expect that for My < +3.5 our density estimates are too 
low. However, we are interested mostly in the luminosity 
functio n for fainter co ntact binaries where the Hipparcos 
sample jRucinskl 2002) ran out of low luminosity objects. 
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Table 2. Table 2: The luminosity function for ASAS sample of 
EW systems with log P < —0.25. n is the number of EW systems 
in a 1-mag. wide bins centred at Mv within the distances ri and 
r-2- ri corresponds to V = 8 and My — 0.5 while T2 corresponds 
to V = 13 and My + 0.5. <j> is the luminosity function in units of 
10~ 7 stars/pc 3 ; its error e<j> is evaluated from y/n. 



M v 


n 


ri (pc) 


>'2 (PC) 


<t> 


£0 


2.0 


489 


199.5 


1258.9 


0.80 


0.04 


3.0 


72<> 


125.9 


794.3 


4.74 


0.18 


4.0 


564 


79.4 


501.2 


14.61 


0.62 


5.0 


229 


50.1 


316.2 


23.62 


1.56 


6.0 


26 


31.6 


199.5 


10.67 


2.09 



ASAS x 3 
HIP 

MS / 500 



1 2 3 4 5 6 7 

M v 

Figure 5. The luminosity function for the EW systems with 
logP < —0.25 in the ASAS survey (the continuous line) is com- 
pared with the results based on the Hipparcos sample extending 
to V m ax < 7.5 (the broken line) and with the solar neighbour- 
hood data scaled down by factor 500 X (the dotted line). Note 
that for the direct comparison with the Hipparcos results, the 
ASAS counts have been multiplied by a factor of 3x. The error 
bars show the mean standard errors resulting from the statistical 
uncertainties only and do not include the systematic uncertainty 
of the ASAS detection correction factor estimated at about 30%. 



The ASAS photometric survey covered the 5-magnitude 
interval 8 < V < 13. The two magnitude limits, expressed 
as distances, are shown as lines in Figure 2] The other lim- 
itation comes from our My calibration (Section [3J which is 
valid within 1.5 < Mv < 5.5. For each l-mag.-wide bin in 
My, we found the near and far limits (ri and r2 in Tabled, 
and simply counted the systems in the respective volume 



4tt/3 x C x (ri 



assuming that the survey South of 



S = +28° covered C = 0.735 of the whole celestial sphere. 
The results are given in Tabled while the luminosity func- 
tion is shown in Figure |S] 



5 THE LUMINOSITY FUNCTION 

Figure compares the ASAS luminosity function, after cor- 
rection by 3 x (for systems missed because of the amplitude 
detection selection and of the removal of systems with P > 
0.562 days), with the Hipparcos luminosity function and 
with the Main Sequence luminosity function iWielen et alJ 
Il983t) scaled down by a factor of 500 times. The ASAS and 
Hipparcos sample results are complementary, with the ASAS 
sample improving the luminosity function within the inter- 



val 3.5 < Mv < 5.5 where the Hipparcos sample had very 
few objects. While the luminosity function is better defined 
in the statistical sense with hundreds of binaries per AMy 
bin, the correction for the systematic effects by a factor of 
(3 ± l)x leaves a rather wide margin for improvement. 

The comparison of the luminosity functions in Figure |S] 
confirms the previously derived relative spatial frequency 
of contact binaries in the solar neighbourhood of 0.2%. Al- 
though the relative frequency of bright contact binaries, 
Mv < +3.5, appears low at about 0.1%, the lower num- 
bers of intrinsically bright binaries may be caused by our 
neglect of the IS reddening and extinction and by the in- 
applicability of the simplified Mv calibration. The relative 
frequency reaches 0.2% for the interval +3.5 < Mv < +5.5, 
but then drops to 0.1% in the Mv = +6 bin, and to zero 
for still fainter objects. The sharp decline in the numbers of 
faint systems is related to the well established period cut- 
off at 0.215 - 0.22 days 3 which still remains unexplained. It 
should be noted that the relative frequency may be higher 
in the central galactic disk at about 1/130 iRucinskilll99ctl 
although this particular estimate suffered from photometric 
blending effects in the direction of the Galactic Bulge which 
may have raised the frequency by a factor of about 2x. 
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of B. Paczyhski, D. Szczygiel, B. Pilecki and G. Pojmahski 
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Support from the Natural Sciences and Engineering 
Council of Canada is acknowledged with gratitude. 



REFERENCES 

Eggen, O. J., 1967, MemRAS, 70, 111 

Hamilton, W.C., 1964, Statistics in Physical Science: Es- 
timation, Hypothesis Testing, and Least Squares, The 
Ronal Press Company, Chapter 4. 

Paczyhski, B., Szczygiel, D., Pilecki, B., Pojmahski, G., 
2006, MNRAS, in press (astro-ph/0601026 ) 

P ribulla, T. fc Ru cinski, S. M., 2006, AJ, in press 
(astro-ph/0601610) 

Rucinski, S. M., 1994, PASP, 106, 462 

Rucinski, S. M., 1998, AJ, 116, 2998 

Rucinski, S. M., 2000, AJ, 120, 319 

Rucinski, S. M., 2001, AJ, 122, 1007 

Rucinski, S. M., 2002, PASP, 114, 1124 

Rucinski, S. M., 2004, New Astr. Rev., 48, 703 

Rucinski, S. M., & Duerbeck, H.W., 1997, PASP, 109, 1340 

Weldrake, D.T.F., Sackett, P.D., Bridges, T.J., & Freeman, 
K.C., 2004, AJ, 128, 736 

Wielen, R., Jahrei/3, H., & Kriiger, R. 1983, in "The Nearby 
Stars and the Stellar Luminosity Function," Proc. IAU 
Colloq. 76, eds. A. G. Davis Philip and A. R. Upgren (L. 
Davis Press, Schenectady, N.Y.), p. 163 



3 The system CC Com delineated the cu t-off at 0.2207 d. for a 
long time, but a new discovery in 47 Tuc bv lWeldrake et aI.H2004l) 
has moved it to a still shorter period of 0.2155 d. Shor t periods 
are e xpected and indeed observed in globular clusters jRucinskil 
2000). 



